Study of the Localization Transition on a Cayley-tree via Spectral Statistics 
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The spectral statistics of a Cayley-tree is numerically studied. The statistics are non-universal 
due to the high ratio of boundary sites. Once the boundary sites are connected to each other in a 
way that preserves the local structure of the tree the universal statistics of the spectra is recovered. 
A clear localization transition is observed as function of on-site disorder strength, with a critical 



disorder Wc = 11.441^;^^ and critical index u 
field value, while the value of Wc is puzzling. 

PACS numbers: 



The properties of the Anderson transition have gener- 
ated much interest since it was first predicted 1]. The 
transition is characterized by a lower critical dimension, 
believed to be equal to two, below which all states are 
localized for any amount of disorder. Above the lower 
critical dimension a transition between extended and lo- 
calized states appears at some critical value of disorder 
(or energy). An upper critical dimension, above which 
the transition may be described by a mean field theory, 
is not well established. 

The Anderson transition is usually characterized by 
two parameters. The critical disorder Wc at the middle 
of the band (where W is the width of the distribution 
from which the on-site energies are drawn in the canon- 
ical Anderson model defined in Eq. (Q) and the critical 
index v. The dependence of Wc and v on the dimension- 
ality d has been the subject of many recent numerical 
studies. For ci = 3,4 the values are well established - 
Wc - 16.5, V ~ 1.5 for d = 3 0, and Wc ~ 35, 
j/^lforc? = 40,|^. The mean field value of the crit- 
ical exponent is equal to 1/2. Assuming that the upper 
critical dimension is equal to infinity, an extrapolation 
equation for v 0.8/(c? — 2) -f .5 was proposed Ver- 
ification of this extrapolation was obtained by studying 
the Anderson transition for bifractal system, where it was 
demonstrated that d should be replaced by the spectral 
dimension ds- Similarly, the critical disorder is also ex- 
trapolated by 16.5(d — 2) (again for bifractals d is 
replaced by dg) (J]. 

Our main goal in this study is to identify the metal- 
insulator transition in a disordered Cayley-tree, and to 
study its properties numerically by spectral statistics. 
Although many studies were performed for Cayley-tree 
structures, to the best of our knowledge, no studies were 
performed using spectral statistics Moreover, it is 
not trivial to extend the above extrapolations for v and 
Wc to the Cayley-tree. From analytical calculations it 
is known that for the Cayley-tree v = 0.5 0,0,0) i-e-j 
d — oo. On the other hand, a mobility edge is predicted 
at some finite disorder, which is hard to conciliate with 
the extrapolation formula for Wc which gives Wc = oo if 
(i = cx) is plugged in, and Wc < if the spectral dimcn- 
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The value of u fits nicely to its mean 




FIG. 1: The distribution P(s) for an L = 7 Cayley-tree for 
different W values, with all eigenvalues taken into account. 
The Wigner distribution (Eq. as well as the Poisson dis- 
tribution @ are indicated in the plot. 



sion of a Cayley-tree ds = 4/3 |lO| is inserted. 

We based our calculations on the usual tight-binding 
Hamiltonian, 



H 



i <i,j> 



(1) 



where the left part of H stands for the disordered on- 
site potential. The on-site energies, are uniformly 
distributed over the range —W/2 < Si < W/2. The 
right part is the hopping element which is set to 1, 
and < i,j > denotes nearest neighbors. Here we 
considered a tree where each site is connected to two 
sites below it. We diagonalize the Hamiltonian ex- 
actly, and obtain N eigenvalues Ei (where N is the 
number of sites in the tree) and eigenvectors ipi. The 
calculations are made for K different realizations - 
K = 4000, 2000, 1000,..., 125, 64 for the correspond- 
ing tree sizes: N = 63, 127, 255, 2047, 4095 or L = 
6, 7, 8, ... , 11, 12 (where L is the number of " generations" 
in the tree). 
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FIG. 2: The distribution P(s) for an L = 10 Cayley-tree in 
which the boundary leaves were connected to each other, thus 
preserving the local Cayley-tree structure while avoiding the 
peculiarities introduced by the boundaries. A clear transition 
from Wigner to Poisson is observed as function of disorder. 



FIG. 3: 7 as function of W for different tree sizes L. The typ- 
ical behavior for finite size transition is seen, where a crossing 
in the size dependence of 7 between the metallic (small values 
of W) and locaUze (large value of W) regime is seen. 



We have calculated the distribution P{s) of adjacent 
level spacings s, where s = {Ei+i — Ei)/< Ei+i — Ei >. 
Typical result are presented in Fig. It can be clearly 
seen that there is almost no change in the distribution as 
function of the disorder once the disorder is above W = 1. 

This unusual behavior of the nearest-neighbor level 
spacing can be attributed to the special form of the 
Cayley-tree. Half of the sites in the tree are boundary 
"leaves"- sites at the boundary of the tree which are not 
connected any further. This peculiar structure of the tree 
is known to lead to unusual behavior such as a jump in 
the participation ratio at the mobility edge ■ -f^^f. 
[1^ this peculiarity is remedied by connecting each of the 
boundary leaves randomly to two other leaf sites. Thus, 
the local structure of the Cayley-tree is preserved, while 
there are no boundary leaves. P{s) for such a tree is de- 
picted in Fig. 12 As expected the distribution is shifting 
from the Wigner surmise distribution (characteristic of 
extended states). 



to a Poisson distribution (localized states), 



(2) 



(3) 



We can recognize the Anderson transition also by not- 
ing that all curves intersects at s ~ 2 and the peak of the 
distribution "climbs" along the Poisson curve for larger 
values of W. The transition point can be established 
more accurately, as shown in Ref. 0], from calculating 



7: 



/-P(.)-/-Pa.) 



(4) 



7 ^ as the distribution tends towards the Wigner dis- 
tribution, while 7 — > 1 if the distribution approaches the 
Poisson distribution. One expects that as the system size 
increases, the finite size corrections will become small re- 
sulting in a distribution closer to Wigner distribution in 
the metallic regime and to Poisson in the localized one. 
At the transition point the distribution should be inde- 
pendent of the system size. Indeed, this is the behavior 
seen in Fig. |31in which 7 decreases with system size for 
small values of W while it increases with size for large 
values of W. All curves seem to cross at a particular 
value of disorder signifying the critical disorder. 

From finite size scaling arguments one expects that 
around the critical disorder 7 will depend on the the dis- 
order and tree size in the following way: 



-f{W,L)^j{Wc,L) + C 



W 



1 



(5) 



where C is a constant. This relation enables us to extract 
the critical disorder Wc = 11.44^^0^ and the critical in- 
dex J/ = 0.5llg'Q4. The scaling of the numerical data 
according to Eq. Q is depicted in Fig. ^ Two branches 
corresponding to the metallic and localized regimes are 
clearly seen. The critical index v fits rather well the mean 
field results mentioned above. On the other hand, the 
critical disorder for the Cayley-tree is lower than in the 
three dimensional case. Thus, while from the extrapola- 
tion equation of v one concludes that the dimensionality 
of the Cayley-tree is infinity (as expected on geometrical 
grounds), from the extrapolation equation of the critical 
disorder one concludes that d = 2.7. This value does not 
correspond neither to the geometric dimensionality nor 
to the spectral dimensionality dg = 4/3. 

It is also interesting to check the behavior of the inverse 
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FIG. 4: The scaling of 7 according to Eq. ^ for different 
tree sizes L. Two brancfies, corresponding to the metallic and 
localized regimes, appear. 



regime / 1. The distribution of I as function of the 
disorder is depicted in Fig. [S] For small values of W the 
distribution is peaked at small values of /, as expected in 
the metallic regime, while for larger values of W the dis- 
tribution is very wide. The transition of the distribution 
between the metallic regime and the localized one is, as 
expected, smooth. 

In conclusion, the spectral statistics of a Cayley-tree 
depends strongly on the boundary condition. When the 
boundary leaves are connected to each other in a way that 
preserves the local structure of the tree, a clear Anderson 
transition is observed. As expected the critical index z/ 
corresponds to its mean field value. On the other hand, 
the critical disorder for which the localization transition 
occurs does not fit into the usual extrapolation formulas 
when either the geometric or the spectral dimension of 
the Cayley-tree are used. 

We acknowledge very useful discussions with A. D. 
Mirlin on the importance of the boundary sites, and sup- 
port from the Israel Academy of Science. 




FIG. 5: The distribution of the IPR for different values of 
disorder W for an L = 10 tree. 



participation ratio (IPR) defined as: 



(6) 



In the metallic regime / ^ ^/N, while in the localized 
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